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We study the optical properties of small gold chains doped with different transition metal (TM) 
atoms (Ni,Rh,Fe) by using the time-dependent density-functional theory (TDDFT) approach. The 
optical absorption spectrum of such systems demonstrates a collective plasmon mode when the 
number of atoms is larger than approximately 10, and this mode splits into two peaks when the 
Au chain is doped with some of the TM (Ni,Fe) atoms. We associate the additional peak with 
a local plasmonic mode which corresponds to the charge oscillations around the potential created 
by the d-orbitals of the impurity atoms. The effect takes place when the number of TM atoms is 
much smaller than the number of Au atoms. This behavior is opposite to the case of larger bulk 
noble-metal-TM clusters (radius >lnm), where the doping with TM atoms does not lead to any 
generation of additional modes, and often leads to a suppression of the main plasmon peak. This 
phenomenon of tuning the optical properties of nanosystems with transition metal atoms can be 
used in many practical applications in the field of nanotechnology. 

PACS numbers: 78.67.Bf, 36.40. Gk, 73.22.-f, 71.15.Mb 



Introduction -The physical properties of bimetallic no- 
ble metal-transition metal (TM) nanostructures is an ac- 
tive area of both experimental and theoretical studies, 
due to a great potential for such systems to be used in 
modern technologies, including catalysis [1, biomedicine 
[2], photophysics [3], and photonics [4 (see also Refs. [SJ 
[6] for a general review) . One of the main reasons for this 
is that many noble metal systems, like bulk and layered 
extended Au and Ag, have important properties in vis- 
ible light frequncy range including plasmon excitations, 
and these properties can be tuned by changing both the 
geometry and chemical composition of the clusters, by 
doping them with TM atoms. Since the TM atoms have 
both extended s- and localized d-states, one may expect, 
similar to the bulk and layered TM systems, interesing 
and unusual interplay of the role of these two types of 
states in determining the properties of the systems, in- 
cluding the optical response. 

In the case of noble-metal-TM systems, the structures 
often have core-shell type, with the TM core. For such 
large particles, the classical Mie's theory can used to 
study the absorption spectrum [7 . As calculations show, 
usually the magnitude of the surface plasmon noble metal 
peak gets suppressed when the number of TM atoms in 
the system increases (see, e.g., Refs. [5j[6]). 

Recently, an interest in studying the optical proper- 
ties of small noble atom chains has arisen, in particular 
due to a possibility to assemble the Au atom chains on 
NiAl(llO) surface within the STM[8 . Since the opti- 
cal properties of such chains, similar to other nanoparti- 
cles, would have many practical applications, the ques- 
tion whether such small systems can demostrate a collec- 
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tive optical response naturally arises in this case. In gen- 
eral, the Mie's theory is not applicable when the chains 
are rather short, therefore one needs to use quantum ap- 
proaches, like TDDFT. As it was shown by Kummel et 
al[9], who applied both TDDFT and quantum fluid dy- 
namics in local current approximation (LCA) to small 
Na atom clusters, one can identify the center-of-mass os- 
cilaltions of the electronic charge with respect to posi- 
tive background as plasmons even in the case of rather 
small clusters, similar to the bulk case. Several studies of 
plasmon modes in noble and alkali metal chains within 
TDDFT approach are reported in literature. In partic- 
ular, Lian et. al.[10 studied the possibility of plasmon 
excitations in Au chains of different length and found that 
the collective plasmon mode in the adsorption spectrum 
arises when the number of atoms is larger than approx- 
imately 10. The possibility of central and end plasmon 
peaks in the case of Na chains with an Ag edge atom was 
shown by Yan et al[TT].They also show that in case of Ag 
atom chains d electrons can lead to decrease intensity 
and energy for transverse mode but does not effect much 
the longitudinal modes [T2] .Another interesting feature is 
observed at surfaces by Pohl et.al.[l3] where they see 
acoustic plasmons due to collective effects of surface and 
bulk electrons which can be the case for future study of 
chains on NiAL substrate where Au chain electrons and 
NiAl bulk electrons can lead to such excitations. 

To our knowledge, no attention has been paid to the 
optical properties of TM-doped noble atom chains. In 
order to get some nontrivial different behavior of the 
bimetallic noble metal-TM nanosystem, the number of 
TM atoms must be rather small. On the other hand, one 
can expect a pronounceable effect from these atoms only 
when the number of host (Au) atoms is not too large, 
when the contribution to the density of states from the 
TM atoms is not negligible. Thus, weakly-doped TM- 
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atom Ail chains would be a good candidate for a system, 
where one can expect non-trivial results in the optical 
properties due to TM atoms. 

In this Letter, we apply a TDDFT approach to study 
the optical properties of 2 to 24-atom Au chains weakly 
doped with different TM atoms (Ni, Rh and Fe). We 
demonstrate that in some cases such a doping leads to 
an unusual response of the system, including generation 
of extra plasmon peaks. 

Computational details .-To study the optical absorp- 
tion spectrum of the chains, we applied a TDDFT 
aproach by using the Gaussian 03 code [14] with B3PW91 
hybrid functional [T5} [T6] and LanL2DZ basis set [17]. 
This choice of the potential and the basis set allowed 
us to reproduce reasonably well the experimental results 
for the lowest excited states of the Au dimer [18] [T^. We 
have considered the case of Au n - m X m chains with the 
length up to n=24 atoms, some (in most cases one) of 
which were substituted by 1 < m < 4 X=Ni,Rh and Fe 
atoms. In the case of pure Au chains, our results are in 
a good agreement with other TDDFT calculations [10 . 
Our main focus was the chains with the interatomic dis- 
tance d = 2. 89 A, observed experimentally in the case 
of Au chains on NiAl(llO) substrate[8]. However, the 
cases with different interatomic distances, including the 
optimized Au-Au and Au-X distances, obtained for the 
Au-Au-Au and Au-X-Au trimer chains, were considered. 
As it follows from our calculations, the optical absorp- 
tion spectrum does not depends significantly on the in- 
teractomic distances even in the case of chains with non- 
equidistant nearest-neighbor atoms. Since the computa- 
tional time for a given number of excited states increases 
dramatically with increasing chain size, we considered 
the states below 5eV, in particular in order to study all 
possible absorption peaks in visible range of frequencies. 
This choice allowed as to study the longitudinal plas- 
mon modes in all cases, though we were not able to con- 
sider possible high-energy (ultra-violet) transverse plas- 
mon modes, found for example in Au chains [10 , which is 
an interesting problem for further studies. Once the en- 
ergy levels Ei and the transition dipole moments (i\D\j) 
were obtained, we calculated the absorption spectrum by 
using the standard expression 

A(lj) = Y,(Ei - E 3 mb\j)\ 2 e-^- E ^y^ (1) 

(in arbitrary units), where Y is the peak broadening. In 
order to make the figures more transparent, in plotted 
A(uj) we neglected the contribution of the dipole mo- 
ments less than 0.1 in atomic units, which produce many 
peaks which are much smaller comparing to the main 
peaks 

Results -In the case of Au clusters (bondlength 2. 89 A), 
we have found that a new (collective) longitudinal plas- 
monic mode appears in the optical absorption spectrum 
when the number of atoms in the even numbered chains is 
of order 10 (Fig. la). The position of the peaks moves to 
the infrared region with the number of atoms increasing 
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FIG. 1: The optical absorption spectrum of the Au n (a) and 
Aun-iNi (b) (2 < n < 24) chains. Here and in Figures below 
the interatomic distance is 2.89A 

and becomes close to the asymptotic value ~ 0.6eV when 
the number of atoms is around 20. Since the magnitude 
of the peak grows with the number of atom increasing, 
this shows that this excitation has a collective form, so 
with increasing number of atoms more electrons are in- 
volved in the collective motion. The redshift of the plas- 
mon energy with number of atoms increasing is due to the 
fact that the energy gaps (bandgaps in the many atom 
case) involved in the dipole excitation are reduced in this 
case, and hence less energy is required for the transition 
to excited states (see also Ref. [10 the case of Au chains 
with even number of atoms was also considered). Inter- 
estingly, we have found that in the case of gold chains 
with odd number of atoms there is a smaller sattelite 
peak with the energy ~ 0.1 below the main peak, which 
may be attributed to an additional oscillations related to 
the unbonded extra charge from the "odd" atom. Since 
the study of the difference between the spectra in the case 
of even and odd number of atoms is beyond the scope of 
this paper, we farther consider only the cases with even 
number of atoms. The results for the optical absorption 
spectrum in the case of Au chains with optimized dimer 
bondlength 2. 55 A are very similar to Fig. la, with a small 
shift of the position of the plasmon peak. 

As the doped case, first we have considered the Au 
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FIG. 2: The optical absorption spectrum of the Au n -iNi, 
n=18, 20, 24, chains. 



chains with length from 2 to 24 doped with one TM atom 
X=Ni, Rh and Fe. We have found that in the case of Ni 
atom there is an additional plasmon peak close in energy 
to the main peak (Figlb). Since this peak does not exists 
in the case of short chains, it also must have a collective 
nature. The facts that this peak appears when one atom 
is doped and its position almost does not depend on the 
position of the Ni atom in the chain (except when it is 
very close to the chain edge) suggests that it must have 
a localized nature (see the next Subsection). We have 
found a very similar result with a localized plasmon mode 
in the case of iron atom, while there is no such peak in the 
case of Rh atom. Since in the Rh atom case the highest 
occupied s-orbital contains one electron, similar to the Au 
atom, while in the Ni and Fe case this orbital is doubly 
occupied, the presence of the extra peak may be closely 
related to the charge redistribution in the doped chains, 
namely much more stronger involving in the bonding of 
the d-orbitals in the case of Ni and Fe atoms (see below). 
To check this idea we have made calculations also in the 
case of Ag impurity atom with singly-occupied s-state, 
and we have found no extra plasmon peak, similar to the 
Rh case. 

A more close look at the plasmonic modes in the Ni and 
Fe cases when the number of atoms is close to optimal 
value of long chains with almost no size effect (chains 
with the length between 18 and 24 atoms) shows that 
at such lengths there is change of the plasmon mode 
weights when the larger oscillation strength of the domi- 
nant higher mode (n=18) moves to the lower one (n=24) 
having two (in the case of Ni) or even three (in the case 
of Fe) close modes with similar strength when n is close 
to 20 (see Fig. 2, where the results for AuNi chains are 
presented). Since n=20 with the Ni atom in the mid- 
dle corresponds to two separated chains with the size 
close to the critical value ~ 10, when the collective ex- 
tended plasmonic mode appears, one may assume that 
at this length there appear strong coherent oscillations of 
two chains which make the lower mode dominant. The 



FIG. 3: The optical absorption spectrum of the 
Au24-nNi n (n = 1, .., 5) chains. 



magnitude of the lower peak increases with farther n in- 
creasing. Thus, the second mode can be attributed to 
the local plasmon, with higher energy energy compar- 
ing to the collective plasmon at n < 20. The reason for 
this is that the attractive positive background potential 
is stronger for the local mode (see Figs. 4 and 5 below). 
In the case of small number of atoms (between 10 and 
20) naturally the collective modes of two subchains (with 
length 5 to 10) have higher energy comparing to strongly 
coupled to the background local charge oscillations (see 
Fig. la). The fact that the optical spectrum in the case 
of long chains almost does not depend on the position 
of the Ni atom, except when it is 2-3 atoms close to the 
edge suggests that the charge distributed over 4-5 Au 
atoms is involved in the local oscillations. Fir this reason 
the position of the local peak depends much less on the 
chain length comparing to the collective plasmon, where 
the energy of the plasmon (spectrum gap) decreases with 
the number of atom increasing (at n less than ~ 20) 

We have found that the plasmonic peaks in the optical 
absorption spectrum disappear when the number of the 
impurity atoms is around 5. The results for the case of 
quasi-equidistant Ni impurities far from the edges in the 
24- atom Au chain are presented in Fig. 3. Indeed, in the 
case of 5 impurities, or in the case of five 3-atom and 
one 4- atom separated Au chains, the chains are too short 
to excite either collective or local plasmon excitation (see 
Fig.l). Similar results were found in the clubbed case. In 
particular, in the case of Au\^Ni2 there are two plasmon 
peaks in both separated and clubbed cases, which are 
separated by distance ~ O.leF, but the center of the 
peak positions is shifted in the clubbed and separated 
cases. This result also suggests that there are two kinds of 
excitations in the system: collective motion in subchains 
and some local charge oscillations near the AuNi border. 

Discussion -In order to get a better understanding of 
possible reasons of the plasmon peak splitting, we ana- 
lyze the Mulliken atomic charge distributions in the case 
of Au and single-doped Au clusters. The results for the 
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FIG. 4: The Mulliken atomic charge distribution in the case 
of AuigX(X = Au, Ni, Rh). The impurity atom is put at site 
10 from the left. 



Au2o, AuigNi and AuigRh are presented in Fig. 4 (the 
results for the AuigFe case is very similar to AuigNi). 
Since the charge distribution on the edges is almost iden- 
tical in all three cases and there is no second plasmon 
peak in the pure Au chain, it means that the edge os- 
cillations are not responsible for the second peak in the 
doped chains. On the other hand, there is very strong dif- 
ference in charge distributions near the impurity atom at 
the center of the chain in all three cases. In the case of Ni 
atom, the impurity ion has a very large positive charge, 
while the nearest Au atoms are negatively charged. One 
can model this situation as a jellium electron model in 
presence of one high potential well and two neighbor val- 
leys, and the rest of the chain can be considered as the 
standard positive background in the plasmon problems 
(Figure 5). The plasmon oscillations can be modeled as 
local mode and a collective mode in the case of atom 
with impurities (Fig.5d and 5e). To understand the dif- 
ference between the results for the Rh and Ni(Fe) atoms, 
one can find from Fig. 4 that the height of the potential 
well and the depth of the valleys are much smaller in the 
Rh case. In this case, when the charge from the right 
subchain close to the valley moves to the right but the 
valley potential is not deep enough to move the charge 
back, there is no reason for the local oscillations. To es- 
timate the necessary potential strength, one can assume 
that the restoring energy effects are important when this 
energy is of order or larger of the kinetic energy in the 
system. The kinetic energy of electrons can be estimated 
as twice nearest-neighbor hopping energy. We approxi- 
mate it by the expression for the Slater-Koster parame- 
ter V SS(T = 1.32ft 2 /md 2 multiplied by 2, where m is the 
electron effective mass, which we put equal to the free 
electron mass m e for simplicity and d — 2. 89 A is the 
inter-atomic distance. This gives Eki n — (Actu- 
ally, this energy can be lower since the effective electron 
mass in the gold chain can be much larger than m e ). 
The potential energy is of order of local plasmon fre- 
quency ujip, which can be estimated from the bulk for- 



FIG. 5: Schematic form of the positive charge background 
potential for the Auio (a )and AuqX (c) chains. In the case 
of Auio, the plasmon excitation corresponds to the charge os- 
cillation along all the chain (b). In the case of AuqX chains, 
there are possible two kinds of plasmonic oscillations: col- 
lective plasmon excitations in each half of the chain (d) and 
local plasmon oscillations on both sides of the impurity atom 
potential wall (e). 



mula uoi p — 'Anne 2 /m, where n is the electron density. 
To find the value of local density, we can use the ratio 
of the plasmon frequencies in the bulk and in the chain, 
which must be equal to the ratio of the square roots of 
the corresponding densities. Since in the bulk the density 
is 5.9 x 10 22 cm -3 and the plasmon frequency - 5.26eV, 
and the chain plasmon frequency is ~ 0.66eV, this gives 
the chain charge density 7.68 x 10 20 cm -3 , or 1 electron 
per atom. If one models the chain as a tube of some ra- 
dius and length N at0 ms x 2. 89 A, this density corresponds 
to the tube radius 12A. In the case of the doped chain, 
when the impurity charge is shifted to the nearest Au 
atoms, the local charge density increases, which leads to 
an increase of uo\ v . From the estimations above, in order 
to get the local plasmon oscillations one need to have 
ujpi ~ Ek in °° leV, which corresponds to the following 
number of electrons per atom: (leV/0.66eV) 2 ~ 2. In 
other words, the nearest to the impurity atoms must get 
one extra electron in order to generate local plasmons. As 
we see from Fig. 4, this analysis gives semi-quantitatively 
correct results in the case of Ni and Rh atoms. In the 
first case, the extra charge is ~ 0.5 electron, i.e. of order 
of one electron, while in the case of Rh atom it it almost 
order of magnitude smaller. 

One can analyze the reason of the charge redistribu- 
tion from the bonding point of view. Since there is only 
one valence s-electron in the case of Rh atom, this elec- 
tron participates in <j bonding, while the other eight d- 
electrons are not very chemically active. Thus, one can 
assume that the charge of the Rh core is +1, similar to 
the Au atoms and the positive and negative charges of 
the Au n -\Rh and Au n "plasmas" are almost the same, 
which results in the same type of the excitations, i.e. one 
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plasmon. In the case of Ni and Fe atoms, the valence 
s-state is doubly occupied which means that mainly d- 
states contribute to the atomic bonding. A more accu- 
rate look shows that d 2 z , d xz and d yz orbitals donate most 
of the charge to the bonding, which results in a positively 
charged core with charge Z > +1. This is the main rea- 
son for the potential well at the impurity atom. Thus, 
one can estimate which sorts of atoms may lead to gener- 
ation of localized plasmon looking for the bonding charge 
redistribution. Obviously the estimation presented above 
is not very accurate, and more refined estimations can be 
easily obtained by improving this one, in particular by us- 
ing correct electron effective mass. However, we believe 
that even such simple estimations give correct order of 
magnitude of the energy of possible plasmon excitations. 

To summarize, we have shown that a weak doping of 
Au chains with some TM atoms may lead to generation 
of local plasmonic modes. This effect is a result of a del- 
icate balance between the s- and d-states of TM atoms, 
which may lead to a strong charge redistribution in the 
chain near the impurity atom. Namely, the double occu- 



pancy of the top s-orbital and partial occupancy of the 
d-states are necessary in this case. We have also proposed 
a simple criteria for possibility of the local plasmon gen- 
eration based on the charge redistribution between the 
impurity and nearest host atoms. These results show 
that some TM atoms can be used to tune the optical 
properties of nanostructures, which can have many pos- 
sible areas of applications: from solar cells to sensors. 
Several important question related this problems remain 
open, including the role of the substrate in the optical re- 
sponse of the chains and possibility of magneto-plasmons 
generated by ferromagnetically-ordered TM impurities 
(see, e.g., Ref.[20 ). We are planning to study these prob- 
lems in nearest future. 
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